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ELIMINATION  OF  SHORT  PULSE  INTERFERENCE  IN 
PULSE  COMPRESSION  RADARS 


INTRODUCTION 

The  Radar  Division  of  NRL  is  currently  studying  environmental  mapping  from  the  viewpoint 
of  improving  automatic  detection,  tracking,  and  radar  control  processes  [1] .  An  important  factor 
in  this  study  is  the  availability  of  data  from  a  coherent  pulse  compression  system.  This  radar  is 
broadband  with  coherent  processing  and  pulse  compression.  The  data  recorded  will  be  /  and 
Q-channel  voltages  after  pulse  compression.  This  means  that  one  environmental  factor  we  cannot 
examine  in  these  data  is  short  pulse  interference  because  the  pulse  compression  filter  spreads  this 
signal  in  time  by  the  uncompressed  pulse  width.  This  report  describes  an  analysis  of  this  phenomena 
by  simulating  the  operation  of  the  signal  processing  hardware  [2]  on  the  computer-generated  radar 
signals.  The  principal  features  of  the  simulation  will  be  (1)  simulating  radar  returns,  (2)  simulating 
a  pulse  compression  network,  (3)  performing  MTI  cancellating,  and  (4)  performing  CFAR  detection. 

Basically,  the  simulation  was  developed  in  the  following  manner.  First,  a  linear  FM  target 
return  was  generated  and  passed  through  a  pulse  compression  network.  The  next  step  was  to  place 
a  Hamming  weighting  function  on  the  filter  and  observe  the  effect  of  this  on  the  sidelobes  of  the 
compressed  pulse.  These  features  were  examined  by  plotting  these  functions.  The  plots  generated 
are  discussed  in  greater  detail  in  a  following  section.  The  next  step  was  to  add  noise  to  the  wave 
form.  The  ratio  between  target  and  noise  was  set  to  yield  a  25  dB  signal  to  noise  (S/N)  ratio  after 
passing  through  the  pulse  compression  network.  Then,  short  pulse  interference  (SPI)  was  added  to 
the  signal.  The  level  of  the  SPI  was  set  40  dB  higher  than  the  average  noise.  The  next  major  develop¬ 
ment  was  a  loop  to  generate  two  sequential  PRFs  (one  with  SPI,  one  without)  which  required 
simulating  the  Doppler  shift  of  the  target  and  the  mechanism  to  detect  the  SPI.  Since  this  was  a 
simulation,  the  software  used  to  detect  the  SPI  can  also  be  used  as  part  of  an  MTI  used  for  detection 
after  pulse  compressions.  Ground  clutter  set  at  the  same  amplitude  as  the  SP.  was  generated  to  test 
the  MTI. 


SIMULATION  OF  TARGETS 

The  waveform  generated  was  one  achievable  by  an  L-band  surveillance  radar  with  pulse  com¬ 
pression.  The  waveform  was  simulated  at  baseband  because  of  the  large  number  of  samples  required 
to  obtain  a  faithful  representation  of  the  signal  at  IF  or  RF. 

The  return  from  an  isolated  point  target  was  generated  from  the  following  expressions  for  the 
in-phase  Xt(t)  and  quadrature-phase  Xq(t)  signals 
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where 

T  =  120  ps,  the  uncompressed  pulse  length, 

H  =  2n  KTt), 
t r  =  3 14159, 

r  =  0.5  ps,  the  compressed  pulse  length, 

=  7r/4,  an  arbitrary  phase  term,  and 
A  =  signal  amplitude. 

These  functions  are  sampled  over  the  uncompressed  pulse  length  of  T  at  points  At  apart  where 
At  =  r/4. 


DESCRIPTION  OF  FILTERS 


The  matched  pulse  compression  filter  (PCF)  for  the  waveform  described  in  Eq.  (1)  is  equiv¬ 
alent  to  a  correlation  receiver  which  correlates  the  received  waveform  with  the  complex  conjugate 
of  the  transmitted  waveform.  This  complex  conjugate  is  given  by 


tpf2  fiTt  \ 

T  2  ) 

(lit2  uTt  \ 

T'~) 

The  effect  of  passing  the  returned  waveform  through  this  filter  is  described  by 

Yj(t)  -  Xj(t)  F ^ j(t)  +  X Q(t)  F^Q(t) 

FQ(t)  =  -XQ(f)  FMI(t)  +  Xj(t)  FMQ(f), 


(2) 


(3) 


where  X*F  is  the  cross  correlation  of  X  and  F,  i.e., 

Y(t)  =  fX(u  -  t)F(u)  du 
u 

and  the  asterisk  represents  the  conjugate.  The  function  Z(f)  =  Yt(t )2  +  Yg(t)2  is  plotted  in  Fig.  1 
in  dB  form  relative  to  the  largest  value. 


APPLICATION  OF  FILTER  WEIGHTING 


To  achieve  better  sidelobe  levels  from  the  PCF,  a  Hamming  weighting  function  was  generated 

from 


W(t)  =  Ah  +  (1  -  A h )  cos 


t  -  772  \ 
T  )' 


(4) 
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where  Ah  =0.54,  and  applied  to  Eq.  (2)  to  obtain 


and 


which  will  be  used  as  the  pulse  compression  filter  from  now  on.  The  result  of  operating  on  the 
signal  of  Eq.  (1)  in  a  manner  analogous  to  Eq.  (3)  is  shown  in  Fig.  2.  While  the  main  beam  is  broad¬ 
ened,  the  first  sidelobe  is  lowered  nearly  30  dB. 


GENERATION  OF  NOISE 

Gaussian  distributed  noise  samples  may  be  generated  from  uniformly  distributed  random  num¬ 
bers  by  the  relationship 


Gj  =  OyJ-2  log  (t/j)  cos  (2kU2)  ,  (5a) 

where  o  is  the  standard  deviation  desired  and  U j  and  U2  are  random  numbers  uniformly  distributed 
on  the  (0,1)  interval.  An  independent  Gaussian  variable  can  be  generated  from  the  same  two  ran¬ 
dom  numbers  by 

Gq  =  0\J-  2  log  (U^)  sin  (2irU2)  .  (5b) 
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Fig.  2  —  Response  of  a  Hamming  weighted  pulse  compression 
filter  to  a  noise-free  target  signal 


Adding  these  Gaussian  noise  terms  to  the  waveform  of  Eq.  (1),  and  setting  A  (of  Eq.  (1))  to 
0.91287  and  o  (of  noise)  to  1.0,  we  generate  a  postpulse  compression  S/N  of  approximately  25  dB. 
This  ratio  was  observed  by  passing  signal  alone  (Eq.  (1))  through  the  filter  and  then  by  passing 
noise  alone  (Eq.  (5))  through  the  filter.  The  waveform  consisting  of  signal  plus  noise  prior  to 
filtering  is  plotted  in  Fig.  3,  and  the  result  of  passing  the  waveform  through  the  filter  is  shown  in 
Fig.  4. 


DESCRIPTION  OF  SHORT  PULSE  INTERFERENCE 

Short  pulse  interference  (SPI)  was  assumed  to  be  of  2  ps  duration  and  arbitrarily  placed  about 
1/3  of  the  way  through  the  120  ps  pulse.  This  pulse  length,  which  is  approximately  equivalent  to  a 
radar  range  extent  of  305  m  (1000  ft),  was  chosen  as  typical  of  long-range  surveillance  radars.  The 
SPI  was  set  40  dB  above  the  noise  level.  The  waveform  of  the  signal  plus  noise  plus  interference  is 
plotted  in  dB  form  in  Fig.  5,  and  the  response  of  the  filter  is  shown  in  Fig.  6. 


DESCRIPTION  OF  DOPPLER 

The  next  step  is  to  generate  sequential  pulses  to  detect  the  presence  of  short  pulse  interference. 
(Unless  the  interfering  source  has  exactly  the  same  pulse  repetition  frequency  or  multiple  thereof, 
the  interference  will  not  occur  at  the  same  range  in  each  pulse.)  At  this  time  it  was  also  decided  to 
include  the  effect  of  target  velocity  on  the  returned  waveform.  The  Doppler  component  of  phase 
shift  of  the  waveform  is  given  by 


Pd(t)  = 


4  ”Vtf0 
c  prf 


(6) 
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Fig.  5  —  Typical  target  signal  plus  noise  plus 
interference  before  pulse  compression 


Fig.  6  —  Response  of  pulse  compression  filter 
to  target  signal  plus  noise  plus  interference 
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where 

Vt  =  475  m/s 

fa  =  1000  MHz,  carrier  frequency 

c  =  300  X  106  m/s,  speed  of  electromagnetic  propagation 
prf  =  333  pps,  pulse  repetition  frequency 

=  portion  of  phase  which  is  time  independent. 

We  may  ignore  i//c  because  it  can  be  absorbed  into  the  general  arbitrary  phase  term  of  Eq.  (1). 


DETECTION  OF  SHORT  PULSE  INTERFERENCE 

To  detect  short  pulse  interference,  two  small  routines  were  generated.  The  first  routine  takes 
the  sample  by  sample  difference  between  two  digitized  functions.  The  second  routine  performs  a  CF  AR 
detection  operation  on  a  digitized  function.  The  CFAR  senses  the  environment  and  sets  a  detection 
threshold  based  on  the  mean  and  standard  deviation  of  the  signal  in  a  region  where  the  test  for  a 
detection  is  taking  place.  Specifically,  the  samples  in  a  range  interval  about  the  sample  to  be  tested 
are  used  to  estimate  the  mean  and  standard  deviation  with  the  samples  immediately  around  the 
test  cells  ignored.  For  detecting  short  pulse  interference,  first  the  /-  and  Q-channel  voltages  of  pulse 
two  are  subtracted  from  the  J-  and  Q-channel  voltages  of  pulse  one  and  then  the  envelope  of  the 
difference  is  formed.  This  coherent  subtraction  is  not  necessary  ;  in  fact,  there  are  valid  reasons  for 
performing  the  process  on  the  amplitudes  of  the  pulses.  However,  this  is  what  was  done,  and  the 
envelope  is  described  by  Eq.  (7). 

E(t)  =ypf/(t,l)-  *7(f,2)J2  +  JxQ(f,l)-  *0(t,2)J2  .  (7) 

This  envelope  is  passed  through  the  detector  with  each  sample  being  tested.  Sixteen  samples  on  each 
side  of  the  test  sample  are  ignored,  and  the  next  32  samples  on  each  side  of  the  test  cell  are  used  for 
reference.  Sixteen  samples  were  blanked  because  that  was  the  expected  maximum  length  of  the 
short  pulse  interference.  If  longer  pulses  were  expected,  then  more  samples  would  have  to  be 
blanked  to  keep  the  signal  we  are  trying  to  detect  from  contaminating  the  reference  cells.  The  trade¬ 
off  is  that  we  wish  to  keep  the  reference  cells  as  representative  of  the  neighborhood  of  the  test  cell 
as  possible,  but  the  reference  samples  must  not  be  contaminated  by  signals  which  extend  in  range. 

The  threshold  is  set  with  the  formula 


T=me+Koe,  (8) 

where  me  is  the  estimated  mean  level  of  the  environment,  and  K  is  the  number  of  estimated  standard 
deviations,  og,  of  the  environment  necessary  to  ensure  a  given  probability  of  false  alarm,  fta.  To 
determine  K,  the  distribution  of  the  environment  was  assumed  to  be  Rayleigh  (as  it  would  be  for 
envelope  detected  independent  noise  samples),  and 

F  X  X 2 

Pfa  =  J  ~  eXP-~  dx  <9> 

t  a2  2cr 
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was  solved  for  T,  where  o  is  the  standard  deviation  (s.d.)  of  the  underlying  Gaussian.  This  yields 
the  equation 


T  =  av/2  1og(l/F/.a)  .  (10) 

The  mean  and  s.d.  of  a  Rayleigh  distribution  are  related  to  the  s.d.  of  the  underlying  Gaussian  by 


2 


a  =  a  /  2 
r  '  2 


Substituting  Eq.  (11)  into  Eq.  (8)  for  the  estimated  values  of  the  Rayleigh,  we  obtain 

T 


(11) 


(12) 


Solving  Eqs.  (12)  and  (10)  simultaneously  we  find  that  a  can  be  factored  out  leaving 

/21og(l/PfQ)-/' 


K  = 


(13) 


'2~2 


For  a  typical  value  of  Pfa  of  10'  6,  we  find  K  =  6.11.  This  K  value  will  actually  yield  a  Pfa  >  10  6 
since  the  calculation  of  K  contains  the  assumption  that  the  estimates  are  the  true  values. 

With  the  large  interference  to  noise  level  quoted  (40  dB),  the  interference  was  easy  to  detect. 

It  is  expected  that  if  the  interference  level  were  low  enough  that  it  could  not  be  detected,  then  after 
the  energy  in  the  short  pulse  was  spread  over  the  entire  range  interval  by  passing  through  the  pulse 
compression  filter,  then  it  would  be  too  small  compared  to  the  range  compressed  target  pulse  to 
cause  any  serious  problem.  Once  the  samples  which  contain  interference  have  been  determined, 
these  samples  are  blanked  in  both  pulses.  Both  pulses  are  blanked  because  we  intend  to  use  pulse 
to  pulse  coherent  processing  later.  The  response  of  the  filter  to  the  signal  waveform  when  inter¬ 
ference  is  present  but  blanked  is  shown  in  Fig.  7.  Close  examination  of  Figs.  7  and  4  (when  inter¬ 
ference  was  not  present)  reveals  only  minor  differences,  indicating  that  blanking  2  jus  of  a  120  ps 
record  has  little  effect  on  the  filter  response.  The  difference  in  dB  levels  in  Figs.  7  and  4  reflects 
only  the  fact  that  the  dB  values  are  relative  to  different  reference  values  in  each  case. 


GENERATION  OF  CLUTTER 

The  clutter  simulated  in  the  model  is  ground  clutter  because  the  relatively  low  PRF  (333  pps) 
of  the  system  cannot  cancel  clutter  with  larger  Doppler  spreads  such  as  rain  clutter.  The  clutter 
return  was  simulated  by  taking  a  large  number  (in  our  case,  20)  of  point  targets  with  random  ranges 
(in  the  vicinity  of  the  true  target),  random  Dopplers  (between  ±0.1  m/s)  and  random  phases.  The 
amplitude  is  set  to  yield  an  average  clutter  value  equal  to  the  short  pulse  interference  value.  The 
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Fig.  7  —  Response  of  pulse  compression  filter  to  target  signal 
plus  noise  plus  interference  with  interference  blanked 


waveform  consisting  of  target  return,  noise,  interference,  and  clutter  for  one  pulse  is  shown  in  Fig.  IS 
The  second  pulse  without  the  interference  is  practically  identical  because  the  noise  and  target  levels 
are  so  far  below  the  clutter,  the  clutter  decorrelates  very  little  pulse  to  pulse,  and  the  interference 
on  the  same  level  as  the  clutter  occupies  such  a  small  portion  of  the  pulse.  The  basic  similarity 
between  the  pulses  is  shown  in  Fig.  9  which  shows  the  envelope  discussed  earlier  in  Eq.  (7),  i.e..  the 
difference  of  pulses  one  and  two.  Comparing  Fig.  9  to  Fig.  5,  we  obtain  an  estimate  of  the  clutter 
residue  of  the  \1TI,  and  there  is  roughly  a  6  dB  difference  in  interference  to  noise  level  in  the  two 
figures.  If  3  c!B  of  the  difference  is  the  result  of  processing  the  noise  through  the  MTI,  then  ihe 
clutter  residue  is  3  dB  above  the  noise.  The  original  clutter  level  was  approximately  40  dB  above 
the  noise. 

The  operation  of  the  filter  on  the  waveform  of  Fig.  8  is  shown  in  Fig.  10.  We  know  from  Fig.  6 
that  the  interference  can  mask  the  signal,  and  now  we  see  that  clutter  even  covers  up  the  distinctive 
interference.  If  we  take  the  range  filter  output  of  both  pulses  and  pass  them  through  the  MTI,  we 
see  that  the  clutter  cancels  and  we  now  observe  the  response  to  the  short  pulse  interference,  Fig.  1 1 . 

We  now  have  all  the  types  of  signals  simulated  that  we  set  out  to  simulate  and  all  the  process¬ 
ing  tools  necessary  to  operate  on  these  signals. 


DETECTION  OF  TARGETS 

There  are  two  detection  processes  we  could  use  on  two  pulses.  The  first  is  a  coincidence 
detector  based  on  two  pulses.  The  other  is  a  coherent  process  using  the  MTI. 

Using  the  coincidence  detector,  we  assumed  a  noise  only  environment  with  no  clutter  and 
no  SPI.  Since  we  assume  we  are  detecting  in  a  noisy  environment,  the  probability  of  false  alarm  is 
independent  pulse  to  pulse  and  we  may  attain  a  Pfa  of  10  (i  by  setting  the  threshold  to  achieve 
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Fig.  8  —  Target  signal  plus  noise  plus  interference 
plus  clutter  before  pulse  compression 
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Fig.  9—  Difference  between  pulse  containing 
target,  noise,  interference  and  clutter  and 
next  pulse  containing  target,  noise,  and  clut¬ 
ter  but  no  interference  before  pulse  com¬ 
pression 


Fig.  10  —  Response  of  pulse  compression  fil¬ 
ter  to  signal  plus  noise  plus  interference  plus 
clutter  (Fig.  8) 


Fig.  1 1  —  Difference  in  pulse  compression  fil¬ 
ter  response  to  a  pulse  including  interference 
and  the  response  to  a  pulse  without  interfer- 
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a  Pfa  of  10"3  on  each  pulse.  This  is  approximated  by  a  value  of  K  =  3.76  in  Eq.  (8).  This  detector 
detects  signals  in  a  noisy  environment  but  is  the  wrong  detector  in  a  clutter  environment.  Since  the 
clutter  changes  very  little  from  pulse  to  pulse,  a  sample  which  exceeds  the  threshold  on  pulse  one 
has  very  high  probability  of  exceeding  the  threshold  on  pulse  two.  The  detections  are  not  indepen¬ 
dent  events.  It  would  be  very  hard  to  determine  the  proper  threshold  without  using  Monte  Carlo 
simulations,  so  we  will  assume  perfect  correlation  of  the  clutter  pulse  to  pulse  which  means  to  ob¬ 
tain  a  coincident  Pfa  of  10-6  we  set  the  threshold  for  each  pulse  to  the  value  of  10'6.  In  this  case 
there  are  no  false  alarms,  but  there  are  no  detections. 

In  the  coherent  detection  process  we  may  first  detect  the  presence  of  SPI  and  eliminate  it  from 
our  uncompressed  pulses.  Each  pulse  is  then  passed  through  the  pulse  compression  filter  to  obtain 
the  results  shown  in  Fig.  12.  This  is  virtually  indistinguishable  from  Fig.  10,  showing  the  dominance 
of  the  clutter  over  the  SPI.  Operating  on  the  two  compressed  pulses  with  the  MTI  yields  the  results 
in  Fig.  13.  The  target  is  easily  detectable  over  the  clutter  residue,  and  the  effect  of  blanking  the 
interference  is  negligible. 


SUMMARY 

We  have  seen  that  short  pulse  interference  (SPI)  has  an  adverse  effect  on  pulse  compression 
networks.  However,  we  have  shown  that  SPI  can  be  eliminated  by  subtracting  two  successive 
pulses  prior  to  pulse  compression.  This  has  the  effect  of  eliminating  clutter,  if  present,  and  high¬ 
lighting  the  SPI  which  would  appear  in  at  most  one  of  the  pulses.  The  SPI  can  then  be  eliminated  by 
blanking  the  cell  in  which  it  appears  and  then  continuing  processing  in  a  normal  manner.  This  would 
include  postpulse  compression  detection  of  targets  either  coherently  or  noncoherently  depending 
on  the  presence  or  absence  of  clutter.  Blanking  the  few  cells  occupied  by  the  SPI  seemed  to  have 


Fig.  12  —  Response  of  pulse  compression  fil¬ 
ter  to  signal  plus  noise  plus  interference  plus 
clutter  with  interference  blanked 


Fig.  13  —  Difference  of  pulse  compression 
filter  responses  to  two  pulses  where  inter¬ 
ference  is  blanked  prior  to  filtering 
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little  effect  on  the  total  information  content  of  the  long  pulse.  In  the  absence  of  SPI,  no  cells  would 
be  blanked  if  the  threshold  of  the  interference  detector  is  set  properly . 

This  seems  to  be  a  promising  means  of  extracting  targets  buried  in  noise,  short  pulse  inter¬ 
ference,  and/or  clutter,  and  we  hope  to  consider  ways  of  implementing  it  in  radars  in  the  future. 
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